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Abstract 

The jet quenching parameter q is analyzed for a quark jet propagating in an 
anisotropic plasma. The momentum anisotropy is calculated at high temperature of 
the underlying quark-gluon plasma, q is explicitly estimated in leading-logarithmic 
approximation by the broadening of the massless quark interacting via gluon ex- 
change. A plasma instability is present. Strong indications are found that q is 
increasing with increasing anisotropy. Possible implications for the saturation scale 
Q s in A — A collisions are pointed out. 
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I. INTRODUCTION 

Jet quenching in Au — Au collisions at RHIC [l|, 0], i.e. the suppression of tt° and 77 
production at large transverse momenta in central collisions when compared with scaled 
measurements in p — p collisions [3] is considered as one of the important observations in 
favour of the production of a plasma state 0, 0]. 

The experimental data from RHIC are used to determine the actual value of the transport 
coefficient q, which controls the radiative energy loss, responsible for jet quenching in the 
induced gluon bremsstrahlung picture @, 0, H, S Ef], Q Q Q H H G3, Q ■ 

Within the kinetic theory framework the coefficient (for an isotropic medium) is estimated 

by 



q = p J d 2 q ± ql 



da 
d 2 q± 



(1) 



where p is the number density of the constituents of the medium, and is the differential 
scattering cross section of the parton (massless quark or gluon) on the medium. 
It is important to keep in mind that q not only determines the energy loss, but it is 
also related to p j-broadening (per unit length) of energetic partons propagating in the 
medium/plasma [fj. 

Because of many theoretical uncertainties in the determination of q its value maybe 
quoted to be in a broad range of 0.5 — 20 GeV 2 J fm, see e. g Il3l. [l9l, [20| . 

It has been found (see e.g. 0, H [H, H Q, 0, H3, UB, H) that the physics of 



anisotropic plasmas differs from that of isotropic ones, because of the presence of plasma 
instabilities in the former. Due to these observations it requires to reanalyze q in the context 
of anisotropic plasmas. Momentum broadening in a homogeneous but locally anisotropic 
high-temperature system for a heavy quark induced by collisions has been discussed recently 



3 II ] (see also [27], ]32|, ]33j). Here the case for massless partons is of relevance. 
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In the following the anisotropic medium is modeled as a high temperature gas. As pointed 



out in 31 



the temperature is kept as a placeholder of the correct hard scale, like the 
saturation scale Q s , present in the anisotropic system. These assumptions allow an analytic 
and tractable treatment, at least for small anisotropies, considered in this paper. In detail 
the main effect is due to the propagation of the dominantly soft gluon in this medium. 



II. SETUP 

In the following momentum broadening (per unit pathlength) for a massless quark in 
the high energy limit (eikonal limit) in a homogeneous but anisotropic high temperature 
plasma is studied. The following kinematical setup is used: the initial energetic (hard) 
quark propagates with momentum p^ along the z-axis, whereas the beam nuclei collide 
along the y-axis, which is the direction of anisotropy, denoted by the three-dimensional unit 
vector n = (0, 1, 0). 

In the limit under consideration momentum broadening q± takes place only in the x — y 
plane, q\ = q 2 x + q 2 



Instead of kinetic theory |32l . I34j q is calculated in the thermal field theory approach [3J 
using the notation and convention of [sjsj], generalizing the expression given by Eq. (CD) to 
include the finite temperature dependence. The coefficient q is expressed by the thermal 
quark self-energy, however weighted by qj_, 

1 

where explicitly 



q = — Tr ^ W (P)] , (2) 



(p) = g 2 c F / 7S4 ?i [V^(PVA> ] • (3) 



d 4 q 

g is the QCD coupling, Cp = (N 2 — 1)/2N C is the colour factor of the quark coupled to 
the gluon. The hard quark with momentum p' = p + q in the medium remains on-shell at 
temperature T = 0, 

£>(p') = 2nfS + (p' 2 ) = 2nf9(p'°)5(p' 2 ) . (4) 
The gluon propagates with momentum g M through the thermal medium, 

A> = -2(1 + f(q )) Im \ v , (5) 

with the Bose-Einstein distribution f(q°) = (exp(g°/T) — l) -1 . denotes the gluon 

propagator, which is treated in detail in Appendix A. 

In the relativistic case and for small momentum transfer, —q 2 /2 = p-p' ~ q±/2, and with 
p ■ q = p°(q° — v ■ q) ~ 0, v = p/p° = (0, 0, 1), one may approximate 

9 2 C F t f d 4 q 



2p° J (2tt 

by 



Im J-^LqI 27r5 + ((j9 + g) 2 )(l + /(g°))Tr ^(j> + , (6) 
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q = -g 2 C F Im 



W) 4 



2 2T 

gU 



27T(5(g -{7-g) 



p ■ A ■ p 



0^2 



(7) 



when q° << T and q°/p° — > 0. Concerning kinematics: v ■ q = \q\cos9 pq = q 



\q\ 2 sm 2 6 pq = \q\ 2 (l 



x 



with x = q /\q\- 



As already pointed out in 3l| this equation ([7]) tells us that all the information about 
the medium, either isotropic or anisotropic, is contained in the imaginary part of the gluon 
propagator (Appendix A). 



III. HARD-THERMAL-LOOP SELF-ENERGY IN AN ANISOTROPIC PLASMA 



In the hard-thermal-loop approximation [37| the retarded gauge-field self-energy is given 
by H 



U^(q)=g 2 J 



d 3 p u dn(p) 
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u/3 



v v qP 



(2tt) 3 ~ dpP q-v + iej ' ^ 

where f M = (l,p/\p\) is a light-like vector describing the propagation of a plasma particle in 
space-time. The self-energy is symmetric, IP^(g) = IP M (g), and transverse, q l Jl^ u {q) = 0. 

In order to determine IP^ the phase space distribution function has to be specified. The 
following form is used, 

f (9) 



n(p) 



p 2 + £(p ■ n) 2 



Thus, n(p) is obtained from an isotropic distribution 7ijgo([p|) by removing particles with a 

In a more general approach the distribution 



large momentum component along n [23|, [24], [26 
n iso need not necessarily be thermal, as assumed here. 

The functions a, /3, 7 and S, introduced in the Appendix A in Eq. ( 1381) are obtained from 



Eq. (jHJ) • Explicit expressions maybe found in |23l . [2J, [31 



A. isotropic case 



For reference the case for £ = is first considered, for which the propagator is given by 
Eq. (1521 . containing the unscreened transverse mode and the Debye screened lon gitu dinal 
one. Therefore dynamical screening due to Landau damping j36|, as e.g. done in 34], has 
to be included. 

To leading logarithm order (LL) in In ^ ~ In - - in the following the only limit of interest 
- one may keep only terms linear in x in the denominator of the propagator in Eq. (152]) , and 
approximate 



Im 



p ■ A ■ p ^u 2 p ■ B ■ p/q 2 



n 



T 



with uj = g°, p^ 1 = p^/p , by 

Im 



x 2 ) 



+ 



q 2 -Tli 



1 



q 2 — ilmUx q 2 + m 2 D — ilmU-L 



(10) 



(11) 
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Im 



1 



+ 



1 



(12) 



-q 2 + iimi 2 D x/4: q 2 + m 2 D + mm 2 D x/2 
using Eq. (1541) . 

Performing the integrations to LL order with g = |g| < T gives finally for a hard quark 

jet 



fee = : m ■ 



(13) 

Z7r mo 

which coincides with the kinetic result Eq. flBTj) given in Appendix B. Dynamical screening 
indeed guaranties this finite answer. 



B. anisotropic case 

In the hot anisotropic system (£ ^ 0) the gluon propagator of the form of Eq. fT49|) is the 
relevant one, such that 



p-A-p = A y 



1 - x 2 



x 2 q 2 y 1 

1 - ah 



+ A 



G 



x 2 (q 2 - a - 7 ) + ~-\{^ 2x 2 q y 6 



, (14) 



with the transverse momentum component q y = q ' ■ n/\q\ = q y /\q\ into the direction of 
anisotropy. 

First, the contribution of to q of Eq. (JTj) is considered, which even to LL accuracy 
shows the possible presence of the plasma instability. The contribution reads 



9 2 CfT f^l-x 2 



1 - x 2 



x 2 q 2 



I(x, a) 



(15) 



noting that x = v ■ q/\q\ does depend only on angles, as does q y , which are integrated by the 
angular integration in f|T5|) . 

In performing the |g| = q— integration 



I(x,a) = Im / q 3 dq Aa(x, q,a) = — Im / dq — 

Jo Jo q l 1 — 



q 



x 2 ) + Rea + ilma 



(16) 



it has to be kept in mind that Ima is non-vanishing because of Landau damping, however, 
Ima oc x for x — > 0, and that Rea is negative for £ > in a range of x values, \x\ < 1, as 
e.g. shown in Fig. la of 24j. Therefore, in this range, including the point x = 0, there is a 
pole present in the integrand of Eq. (|T6l) . characteristic of the anisotropy. 
Evaluation gives 



I(x, a) 



Ima 



+ 



2(1 -x 2 ) 2 12 
Rea 



1 [T 2 (l - x 2 ) + Reef + (Ima) 2 



a 



Ima 



Ima 



Ima . . 

arctan —— — arctan — h 7rc)(— Rea) 

T 2 (l — x 2 ) + Rea Rea 



(17) 



5 



In the LL approximation this gives 

T , , Ima 
l[x,a) 



1-x 



2\2 



1 , T 7r Rea _ , % 
In h -7 e(-Rea) 

2 m D 2 Ima v ; 



(18) 



where the second term has to be kept, because it reflects a singularity for Ima oc x — > due 
to the anisotropy £ > 0. 

To be explicit, the small £ behaviour is considered, where [24| 



„2 n2 



a= U T (x)+^j(5q 2 -l)m 2 D -^-m 2 D 
1 



+ -n T (x)[3^-l-x 2 (5^-l)]^ , (19) 



which shows that for x — > 

1 



Rea ~ ~Z%ml , (20) 
indeed negative for £ > 0, and 

Ima ~ -\ x(l - x 2 )m 2 D |l + |[3gJ - 1 - x 2 (5q 2 y - 1)]| . (21) 

Following [3l[ the contribution from the first term in Eq. (I18p to is denoted as regular. 
After performing the angular integrations in Eq. (fl5|) it leads at LL order with T » rriu 
to 

& e9 = g^f^ ln— (1 + 0(0), (22) 

with no contribution at first order in £. 

Turning to the regular contribution due to A G one first obtains 



= 9 gJ 3 T I Ml, -rCL - -'• 



„2\ 



, 22 (l-x 2 ) q 2 + a + 1 -^( q 2 -P) + 2q y ~S 

x Im / dg 2 g 2 = -^-^ = . 23 

J x 2 (l-x 2 )(/?-g 2 )g 2 + x 2 (a + 7 )(/?-g 2 )-(5 2 (l-g2) 

The regular LL part reads after the q— integration 

= / dSl q [(l-x 2 )Im^ + 2(l-x 2 )g,Im5+^ilm(a + 7 )] , 

(24) 

with (|2U) and 

Im/3 = -^m^ {1 + £ [2g 2 - 1 - x 2 (3g 2 - 1)]} , 



•^(l-^m?, {1 + | 



Im(a + 7 ) = - T x(l - x 2 )m^ {1 + i[5gj - 3 - a; 2 (7g 2 - 3)]} 



Im5 = ~x(l - x 2 )m 2 D £ (1 - Ax 2 )q y , (25) 
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taken from 2J, |3l| • 



The angular integration in Eq. (1241) then gives 

= 3g 2 C F m 2 D T ^ T_ + 2 
Stt mo 

Summing the two terms Eqs. ( 1221) and ( |26i) the LL transport coefficient in the limit of 
small £ up to becomes 

£Lo = € 9 + ft 9 = 92CFmlT In — (1 + 0(f))- (27) 
For £ = the value for g iso , Eq. ( fTBl . is recovered. This result for g^jo in Eq. (l27|) does 



agree with the result given for k^_ 9 + K r z e9 in Eq. (14) of 31[ when the velocity v is taken 
to be v = 1 for a massless quark jet. Because of this agreement, which is not unexpected, 
one may try to extrapolate %^i so for v — > 1 from the numerical values of k^ 9 z summarized 
in Table 1 of 3l[ even for large values of £ by 



Qaniso ^ ^ ^8) 



almost independent of £ < 100. 



Next the anomalous contribution |3J| due to the second term of Eq. ffTBl is evaluated. In 



LL order only the behaviour for x — > is relevant. With Eq. (1201) it gives 

g 2 C F m 2 D T^ r q 2 



24vr 2 



< J dn,| , (29) 



inducing a logarithmic singularity with x = cos6 pq . The contribution to g£f om starts at 
A short way to see the anomalous contribution is to sum the singular parts in the gluon 



in 



propagator for the static case in terms of the masses, here explicitly quoted from [24 
the limit of small anisotropy £ — > scaled with respect to the Debye mass Eq. (|55|) . i.e. 

rh = m/mD, 

m 2 a = -|(l + cos20„) , 
o 

rh* = | sin 2 9 n , (30) 

and 



m. 



m_ 



l + ^(3cos2# n -l) 
o 



cos 29 „ 



(31) 



where the angle is given by cos# n = q y . 

Since m 2 a and m 2 _ are negative the pole contributions (with the Feynman prescription) 
read 



9 2 C F T 
8vr 3 



Im / dVL n 



q 2 dq 2 



x 



1 



2 2 2 "1 

ml — — 



q 2 + m 2 — ze m\ — m 2 _ q 2 + m 2 _ — ie 



(32) 
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The term of 0(£) comes from the first pole term and leads to Eq. (129|) . whereas the second 
term is obviously of 0(£ 2 ). 

The question about the treatment of the logarithmic singularity in Eq. (|29|) arises. At 
least three possibilities to cut-off the singularity may be discussed: 

(i) One may follow the detailed and plausible arguments given in [3l[ that this soft 
singularity is screened by 0(g 3 ) terms in the gluon propagator, i.e. beyond the HTL ap- 
proximation under discussion. It leads to the replacement of Ima in the second term in the 
denominator of Eq. ffT8l) by Ima ~ x — > x + eg, i.e. it is suggestive to cut the singularity in 
d29D by 

e /^^ 2 £/ -^~2£hi-~2£ln— . (33) 
J x Jo x + eg g run 

This way a finite result is obtained, 

qT m ^ 92CFmlT ^— (34) 
HA 2tt m D 6 ' K J 

which shows a positive, but weak dependence on £ as a sign of the anisotropy for £ > 0.. 

(ii) The origin of the 1/x singularity is traced back to the Bose-Einstein distribution 
f(q°) ~ T/q° in Eq. (jSJ). Pragmatically, in the anisotropic case, this behaviour could be 

modified by q° — > \J (q ) 2 + £(n • <f) 2 = g ^/x 2 + £g 2 , g° > 0. On mass-shell this replacement 
gives the distribution in Eq. (jUJ), and leads to 



--2£ln± (35) 
a; £ 



(iii) To form the anisotropic configuration in momentum space a characteristic time scale 
is present of the order r c ~ 0(l/g£T), for not to large £. It is then natural to cut the 
energies of the constituents in the heat bath by \q°\ > l/r c > g£T, and 

«/7^'4' (36) 

In summary all three options Eqs. (1331 - [36]) lead to a positive contribution of 0(£) at LL 
order to g^o- 



IV. CONCLUSION 

Because of the approximations, keeping only leading logarithmic order terms, a detailed 
quantitative study is not aimed. Qualitatively, however, the main result obtained strongly in- 
dicates that the instability due to the anisotropy leads to q an iso > Qiso- For small anisotropies 
it turns out to be an effect of 0(£), which should be taken into account in the future 
phenomenological comparison with experimental jet quenching data, although the actual 
numerical value depends on the full treatment including the non-leading terms. 

Independent information on q an i SO maybe found in the following references: 



i) A result q an iso > Qiso is found in a numerical simulation in [33j treating jet broaden- 
ing in an unstable non-Abelian (SU(2)) plasma. At early evolution times q an iso — Qiso — 
2.2 GeV 2 / fm, whereas at times when the instability is acting there are numerical indications 
that qaniso increases when compared to cji S0 . 



8 



ii) The transport coefficient q of a fast parton (gluon) propagating in an expanding turbulent 
quark- gluon plasma is estimated in [El, [39, 40]. The coefficient called anomalous in the quoted 
references (but not to be confused with the one in Eq. ({15]) ) is parametrically given by 

q A ~ g^m 2 D T , (37) 

3 

therefore qA/qiso — & 1 9 > 1, even significantly larger than one for small coupling g (assum- 
ing £ not very small). 

These results together indeed show that q an iso > Qiso- 
In A — A collisions the parton saturation scale Q s |41[ is for a large nucleus A related to the 



jet quenching parameter (of gluons) q by Q 2 = 2gV R 2 — b 2 [42, |43_|, |4J, |45J, where R is the 
radius of the nucleus and b is the impact parameter of the collision. Because q is affected by 
anisotropies in A — A collisions, it is not obvious that Q s determined from small— x processes 
in proton/deuteron - nucleus collisions is indeed the same as in A — A collisions. 
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V. APPENDIX 

A. Propagator in covariant gauge in an anisotropic plasma 

Here we summarize the results for the gluon self-energy IP" and the gluon propagator 
zA^ (diagonal in colour) in covariant gauge in an anisotropic medium. 
The self-energy is decomposed as 

IP" = aA^ + fiB 1 ™ + 7C" 1 " + SD^ . (38) 



The tensor basis of IP" is determined explicitly for anisotropic systems in [24j . We use 
here a basis appropriate for the covariant gauge, and define 



,„ K»K U mj q»q u K»K V 

vF = -P^ + — = — g + -Hh + — 39 
K 2 y q 2 K 2 



a 2 K^K V 

B*" = =^=— (40) 

(q-u) 2 K 2 V ' 



= — (41) 



n 



D iw = + K^fii* ; (42) 

where is the heat-bath vector, which in the local rest frame is given by w M = (1, 0, 0, 0), 
and 

q 2 

K* = g M - vP . (43) 

{q-u) 
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The direction of anisotropy in momentum space is given by the vector 

n M = (0, n) . 

Defining the vector 



g M = - (q ■ u)u^ 



and the projector 



introduces the vector 



(44) 
(45) 
(46) 
(47) 

(48) 



where A is the gauge fixing parameter. Upon inversion [46j, the propagator is finally written 
as, with uj = (q ■ u), 



n» = p^n v . 



Both vectors and are orthogonal to g M , and K ■ n — 0. 
The inverse propagator in covariant gauge is then 



1 



a- 1 ) (?, = -<zV + <tf - n^(g, - ^5* 



A"" = A A [A^ - C* v \ + A G 

where 
and 



4 2 

{q 2 - a - l)^jB^ + {uo 2 - (3)C» V + 5^rD^ 



A-/ = (q 2 - a) 



A 5 1 = (g 2 -a- 7 )(a; 2 -^)+5V 2 n 2 



For £ = 0, we recover the isotropic propagator in covariant gauge [37 



= — 



where IT^ = ^IIl and 361 

n T (w,g) 



m 2 } uj 2 



2 f 



(g 2 -n L )g 2 



1 ~~ o R log FT 

2u\q\ uo — \q\ 



(49) 
(50) 

(51) 
(52) 



U L (uo,q) = m 2 D 
In the soft uo limit we remind that 



L) u+\q\ 



log 



2\q\ *uo-\q\ 



- 1 



2 • W 

— m n — tTT——m" J 



(53) 



(54) 



4|g] |g] 2 2|g] 

mn is the isotropic Debye mass, which for a thermal system at temperature T with N t 
colours and Nf light quark flavours is given by 

g 2 T 2 



m 



D 



N f 



(55) 
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Eq. fj4*9]) agrees with the expression given in 43], when the following relations are used 



Q 2 

K» = --^—m» , (56) 
[q-u) 



and 

(q ■ n) ~ 



n» = fit— - 



, (57) 



where the tilde vector in [43] is given by 



<gs = n»~ . (58) 

B. Kinetic theory approach 

In the case of an isotropic medium the transport coefficient q maybe calculated using 
Eq. ([1]), where the differential cross section is related to the spin and colour summed (aver- 
aged) squared matrix element by 

d(J 1 '-Ml 2 , (59) 



d 2 q± 167r 2 s 2 

with s the invariant energy squared. 

When properly taking quantum statistics (Pauli-blocking and Bose-Einstein enhance- 
ment) into account, it is appropriate to modify (PQ) as 

p\Mf^J -0^ [4AT c Ar /n ^(l - ^)|A^|^, + 2(iV^ - l)n BS (l + ^)]7W]y , (60) 

together with the scattering of the quark jet with the thermalised quarks and gluons in the 
hot medium, distributed according to the functions n FD and n BE , respectively. 

At leading order, after screening the effective forward scattering amplitudes j48| by the 
Debye mass l/q±—* l/(q± + mji), the result is 

g 2 C F m 2 D T T 

Q\ kinetic — In , (61) 

2tt rriD 

in agreement with Eq. (|T3l) . This result agrees with the LL one for d < p 2 j_ > jdt = 2 kt 
calculated in [32| (Eq. B32) taking the quark velocity to be v — 1. 
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